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ABSTRACT 

We report the results of a successful 24 hour 6.7 GHz VLBI experiment using the 30 meter radio 
telescope WARK30M near Warkworth, New Zealand, recently converted from a radio telecommunica¬ 
tions antenna, and two radio telescopes located in Australia: Hobart 26-m and Ceduna 30-m. The 
geocentric position of WARK30M is determined with a 100 mm uncertainty for the vertical component 
and 10 mm for the horizontal components. We report correlated flux densities at 6.7 GHz of 175 radio 
sources associated with Fermi y-ray sources. A parsec scale emission from the radio source 1031-837 
is detected, and its association with the y-ray object 2EGL J1032.9-8401 is established with a high 
likelihood ratio. We conclude that the new Pacific area radio telescope WARK30M is ready to operate 
for scientific projects. 

Subject headings: Astronomical Instrumentation 
1. INTRODUCTION 

The New Zealand 30-m radio telescope facility is lo¬ 
cated near the township of Warkworth, New Zealand, 
just 5 km from the Pacific coast. It was built in 1984 as 
a telecommunication antenna (the Earth Station) to be 
used by the New Zealand Post Office, later by Telecom 
New Zealand. In 2010 it was transferred to the Institute 
for Radio Astronomy and Space Research (IRASR) of 
Auckland University of Technology for conversion to a 
radio telescope. 

The 30-m antenna has a special wheel-on-track beam- 
waveguide (BWG) design. This design is often used for 
deep space network and telecommunications (Imbriale 
2005). The feed horn and the receiver are located in 
the pedestal and do not rotate when the antenna moves 
along the azimuth and elevation axes. A system of six 
mirrors directs the beam into the receiver (see Eigure 1). 

The design of the WARK30M antenna should be taken into 
account for data reduction when dual-polarized data are 
processed. 

Conversion of the 30-m antenna took three years and 
required a number of upgrades and repairs. Among them 
was enabling the telescope to track sources with pointing 
accuracy 1', connecting the telescope with a 10 Gbps In¬ 
ternet, installation of an uncooled C-band receiver that 
was donated to IRASR by Jodrell Bank Observatory, in¬ 
stallation of a new Digital Base Band Converter (DBBC- 
2) (Tuccari et al. 2010), installation of a Mark5B+ VLBI 
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recorder, connection to the observatories existing Hydro¬ 
gen maser via a Symmetricom fibre-optic distribution 
amplifier. Detailed description of the upgrades and the 
specification of the Warkworth 30-m radio telescope can 
by found in Woodburn et al. (2015). This work was com¬ 
pleted in late 2014, fringe checks were performed and 
the antenna was ready for its first scientific VLBI exper¬ 
iment. 

As a part of commissioning new antennas, the po¬ 
sition of the antenna reference point should be deter¬ 
mined. A reference point is defined as the point of the 
projection of the movable elevation axis onto the fixed 
azimuthal axis. Eor the analysis of VLBI source imaging 
experiments made in a phase-referencing mode, position 
of VLBI stations should be known with accuracy 10- 
20 centimeters to avoid noticeable image smearing^ (see 
e.g. Chariot et al. 2002)). Eor astrometry applications, 
angular position accuracies of tens of /ias are required, 
making the requirements on the accuracy of station po¬ 
sitions much more stringent: 5-10 mm. One way of es¬ 
timating the position of the antenna reference points is 
through analysis of a combination of a ground survey of 
markers attached to the antennas from a local network 
around the station and GPS observations from the points 
in the local network (see e.g., Sarti, Sillard & Vittuari 
2004; Sarti, Abbondanza & Vittuari 2009). At the time 
of writing this had been planned for the Warkworth 30 m 

^ See http://www.evlbi.org/user_guide/phrefmp.html for an 
example of the effect of position errors in quality of image restora¬ 
tion 
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Beam Waveguide system-view showing 2.5m offset 
between dish and feed axes 



Beam Waveguide system -orthogonal view 
with detail of waveguide mirrors 



Fig. 1.— wark30m beam wave-guide diagram in two orthogonal projections (not to scale). Radio waves from the sky are reflected from 
the main mirror ^1, then reflected from the secondary mirror and guided to the diagonal mirror ^3, which has a circular shape in the 
plane perpendicular to the elevation axis. The parallel beam formed by this mirror travels 2.3 m towards the elliptical mirror ^4 that 
reflects the beam into the waveguide that is tilted at 14.74° with respect to the vertical. The elliptical mirror ^5 sends the beam to the 
diagonal mirror ^6, which directs it into the feed horn mounted in the antenna pedestal. 


antenna but not carried out. An alternative procedure 
for estimation of station position is to use the VLBI tech¬ 
nique to determine group delays and then derive the ref¬ 
erence point position from these group delay measure¬ 
ments. The advantage of the latter approach is that it 
also provides useful diagnostics of the VLBI equipment 
and can be combined with astronomical tasks. 

Here we report results from the first scientific experi¬ 
ment conducted with WARK30M and two Australian radio 
telescopes operated by the University of Tasmania in Ho¬ 
bart and Ceduna (McCulloch et al. 2005) on December 
11, 2014. The VLBI experiment settings are presented in 
section 2. The data analysis for both geodetic and astro¬ 
nomical purposes is described in section 3. Concluding 
remarks are made in section 4. 

2. VLBI EXPERIMENT 

Since the WARK30M C-band receiver has a relatively 
narrow band, approximately 300 MHz, our choice of fre¬ 
quency range for VLBI observations was rather limited. 
We recorded right-circular polarization in a range [6.592, 
6.848] GHz. The signal was split into 16 intermediate fre¬ 
quencies (IFs) of 16 MHz each and recorded with 2 bit 
sampling and an aggregate bit rate 1 Gbps. This setup 
allows us to determine group delay with formal uncer¬ 
tainty 215 ps when the signal to noise ratio (SNR) is 10. 
Two radiotelescopes in Australia participated: ceduna 
and HOBART26 (Figure 2). 

We set three goals for the first scientific experiment: 

• to determine the position of WARK30M with cen¬ 
timeter level of accuracy; 
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Fig. 2.— VLBI network of this experiment. The map shows lo¬ 
cations of the three radio telescopes in New Zealand and Australia 
and the corresponding baselines. 

• to test the effectiveness of the southern VLBI array 
with WARK30M for monitoring radio flux density of 
7 -ray sources; 


• to test the usability of the southern VLBI array 
with WARK30M for association of y-ray sources dis¬ 
covered with Fermi by detecting their parsec-scale 
radio emission. 


The pool of 483 targets included all sources from the 
Radio Fundamental Catalogue^ (Petrov and Kovalev, 

^ The catalogue is available at http://astrogeo.org/rfc 
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2015, in preparation) that are known to be brighter than 
200 mJy at 8 GHz on baselines longer than 5000 km 
and are associated with 7 -ray sources. At the moment, 
20 % of sources with parsec-scale emission brighter than 
200 mJy are known to be 7 -ray loud. Restricting the 
source list only to these objects did not affect our ability 
to measure station positions. In total, 244 target sources 
were scheduled. This includes 7 radio sources from a 
dedicated Australia Telescope Compact Array (ATCA) 
5.5/9.0 GHz survey (Schinzel et al. 2015) never before 
observed with VLBI. These sources are located within 
the error ellipse of Fermi objects but marked as “unas¬ 
sociated” in the the 2 FGL catalogue (Nolan et al. 2012 ). 
As shown by Schinzel et al. (2015), detection of parsec- 
scale radio emission is a powerful method for association 
of 7 -ray sources. Since the number of compact radio 
sources is limited, as a rule of thumb, the probability of 
finding a compact background source with flux density 
from parsec scales brighter than 10 mJy at 8 GHz within 
the error ellipse of a 7 -ray object is less than 10 %. 

Observations were scheduled with software sur_sked^ 
in a sequence that minimizes slewing time and complies 
with the constraints that (a) the next source is located 
at a distance not less than 20 ° and (b) the minimum 
temporal separation between consecutive observations of 
a target source is 3 hours. Target sources were observed 
for 90 s each, except for those seven sources in the vicin¬ 
ity of unassociated 2 FGL objects: they were observed 
for 300 s each. Every hour a sequence of 4 scans of at¬ 
mosphere calibrators was inserted in the schedule: two 
scans of sources at elevation angles in a range [15°,30°] 
and two sources at elevation angles in a range [60°, 88 °]. 
The purpose of including these scans in the schedule was 
to improve the robustness of estimation of residual at¬ 
mospheric path delay in the zenith direction. 

3. DATA ANALYSIS 

The data were transferred electronically to Hobart and 
correlated with the DiFX software (Deller et al. 2007). 
Further processing was undertaken with the software 
package VTMA We determined the complex band¬ 
pass using several strong sources. After applying the 
bandpass calibration, we performed fringe fitting for each 
observation, i.e. we determined group delay, phase delay 
rate, and group delay rate that maximizes the coherent 
sum of visibilities over time and frequency (Petrov et al. 
2011a). The fringe plots for the HOBART26/WARK30M 
baseline were as expected with no obvious anomalies 
(See an example in Figure 3). Unfortunately, the re¬ 
ceiver frequency response at CEDUNA did not match the 
response from other telescope receivers. Fringe ampli¬ 
tude for baselines to CEDUNA at frequencies higher than 
6.620 GHz dropped by a factor of 10 (Figure 4). 

The SNR for baselines to CEDUNA was typically a fac¬ 
tor 4 worse than expected, and, therefore, uncertainty 
in group delays were higher by the same factor. This 
resulted in the non-detections of 19% of target sources 
due to the corresponding sensitivity degradation. Nev¬ 
ertheless, the experiment did not fail, although position 
accuracy achieved was lower than expected. 

^ See http://astrogeo.org/sur_sked 

^ See http://astrogeo.org/web 


3.1. Geodetic Analysis 

Further processing was made with VLBI analysis soft¬ 
ware VTD/Post-Solve^. After evaluation of group de¬ 
lays from visibility data, the theoretical path delays were 
computed using the state-of-the art model and small dif¬ 
ferences between them and the measured group delays 
were formed. A detailed description of this step can be 
found in Petrov et al. (2011b). 

The a priori slant path delays in the neutral atmo¬ 
sphere in the direction of observed sources were com¬ 
puted through numerical integration of differential equa¬ 
tions representing the wave propagation through the 
heterogeneous media. The four-dimensional field of 
the refractivity index distribution was computed using 
the atmospheric pressure, air temperature, and specific 
humidity taken from the output of the Modern Era 
Retrospective-Analysis for Research and Applications 
(MERRA) (Reinecker et al. 2011). The output of the 
model presents the atmospheric parameters at a grid 
0.67° X 0.50° X 72 pressure levels x 6 ^ . A priori slant 
path delays in the ionosphere were computed using global 
ionospheric maps of the total electron content derived 
from Global Navigation Satellite System observations by 
the analysis center GODE. 

In the initial least square (LSQ) solution, positions of 
all stations, except HOBART26, were estimated, as well 
as the coefficients of the expansion for the clock func¬ 
tion and the residual atmosphere path delay in zenith 
direction into the B-spline basis of the 1st degree. Dur¬ 
ing the initial phase of the data analysis, outliers were 
eliminated and the baseline-dependent corrections to the 
a priori weights defined to be reciprocal to formal uncer¬ 
tainties of group delays were determined in such a way 
that the ratio of the weighted sum of squares of residu¬ 
als to their mathematical expectation was close to unity. 
The weighted root mean squares (wrms) of postfit resid¬ 
uals was 114 ps on HOBART26/WARK30M baseline and 
440 ps on baselines to CEDUNA. The increase of resid¬ 
uals with CEDUNA was caused by the mismatch of the 
receiver’s frequency response. 

The final LSQ solution used all VLBI group delays col¬ 
lected from January 1984 through January 2015, in total, 
11.4 million values, including 705 group delays from this 
experiment. Positions of all stations, coordinates of all 
sources, the Earth orientation parameters, as well as over 
1 million nuisance parameters were estimated in a single 
LSQ run. The differences in velocity of WARK30M and 
WARK12M was constrained to zero. Station WARK30M 
has a significant axis offset, the distance between the 
azimuth fixed axis and the moving elevation axis. We 
have adjusted that parameter in our VLBI data analy¬ 
sis. Minimal constraints were imposed to require that the 
net-translation and net-rotation over new position esti¬ 
mates of 44 stations with long history with respect to po¬ 
sitions of observations to these stations in the ITRE2000 
catalogue (Altamimi, Sillard & Boucher 2002 ) be zero. 
This ensures that the positions of all stations, including 
WARK30M be consistent with the ITRE2000 catalogue. 
More details about the parameter estimation technique 
can be found in Petrov et al. (2009). 

The position estimate of WARK30M at epoch 2000.01.01 

^ See http://astrogeo.org/vtd 
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Fringe phase at HOBART26/WARK30M, source 0605-085 



Fig. 3.— Typical fringe phase (Left) and fringe amplitude {Right) at 

Fringe amplutude at CEDUNA/WARK30M, source 0723-008 
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Fig. 4.— Fringe amplitude geduna/wark30m baseline without 
applying the complex bandpass. Receiver bands at these two sta¬ 
tions do not match. 


TABLE 1 

Coordinates and velogities of wark30m derived from 
ANALYSIS OF THE VLBI EXPERIMENT ON EPOGH 2000.0. 


X 

-5115425.60 ± 0.08 

m 

Y 

477880.31 ± 0.02 

m 

Z 

-3767042.81 ± 0.06 

m 

Axis offset 

2.61 ± 0.06 

m 

A 

-19.6 ± 1.2 

mm/yr 

Y 

-2.5 ± 0.5 

mm/yr 

z 

34.1 ± 1.0 

mm/yr 


Since velocity of wark30m was strongly constrained to wark12m, 
the reported velocity is based primarily on geodetic observations 
of WARE 12m station. 

is provided in Table 1. The errors reported in the table 
are the formal uncertainties from the LSQ solution com¬ 
puted in accordance with the error propagation law. Pro¬ 
jected to the local coordinate system, the uncertainties 
are 105, 14, and 12 mm for the Up, East, and North 
components respectively. High uncertainty in the Up 
component is due the necessity to estimate antenna axis 
offset. These two parameters correlate. If we choose not 
to estimate axis offset, but keep its value to 2500 mm, 
as specified by the antenna manufacturer, the position 
uncertainties drop to 38, 12, and 9 mm respectively. 

The largest source of systematic error is the path de¬ 
lay in the ionosphere. In order to evaluate the contri¬ 


Fringe amplitude at HOBART26/WARK30M, source 0605-085 



hobart26/wark30m baseline before applying the complex bandpass. 

bution of the residual ionosphere to site position esti¬ 
mates, we re-processed the CONT-14 VLBI campaign 
with participation of HOBART26 and WARK12M. Dur¬ 
ing 14 days in May 2014 a 17-station array observed a 
geodetic schedule. We processed fourteen 24-hour ex¬ 
periments independently and estimated baseline lengths 
between HOBART26 and WARK12M. We undertook two 
separate analyses. We used ionosphere-free linear com¬ 
binations of group delays at 2.3 and 8.4 GHz in the first 
run. The so-called baseline length repeatability, i.e. the 
wrms of baseline lengths, was 2.7 mm. In the second run 
we used group delays at 8.4 GHz and applied the a pri¬ 
ori ionosphere contribution computed from the GODE 
global ionosphere model using the same algorithm as we 
used for analysis of the experiment with WARK30M. The 
wrms was 6.2 mm, i.e. the contribution of residual iono¬ 
sphere to baseline lengths was 5.5 mm in quadrature. 
Eor deriving the contribution of the ionosphere to the 
vertical component of the WARK30M position estimate, 
we should scale the repeatability of the 2415 km long 
baseline HOBART26/WARK12M by two factors: the geo¬ 
metric factor, the ratio of the Earth’s diameter and the 
baseline length, and the frequency factor: the square of 
the effective frequency observed in GONT-14 campaign 
to the square of the effective frequency observed in our 
experiment. This gives us an estimate of the systematic 
error in the vertical component of WARK30M position: 
50 mm. The systematic error in the horizontal compo¬ 
nent of position vector is smaller by the geometric factor: 
10 mm. This analysis demonstrates that the systematic 
errors due to the mismodeled ionosphere contribution do 
not dominate over random errors. 

3.2. Astronomy Analysis 

Among seven target sources in the vicinity of Fermi 
objects observed in our experiments, one radio source, 
1031-837, has been detected. Its coordinates, determined 
from the global LSQ solution, are a = 10^30’^15^.286 ± 
0L003,(5 = -84°03'08".652 ± 0".004 for the J2000.0 
epoch. No cos (5 factor was applied to the reported un¬ 
certainty in right ascension. 

At the time of our experiment, concurrent measure¬ 
ments of system temperature were not implemented. The 
system temperature at different elevations was measured 
eight days before the experiment, and we found that 
it fits well to the following regression model: Tgys = 
82.6 + 4.3/ sinE, where E is the elevation angle and the 
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temperature is in Kelvins. A priori measurements of the 
system equivalent flux density (SEFD) resulted in the 
value of ^600 Jy. 

Since imaging the data of several scans from a three- 
element array without concurrent Tgys measurements at 
one of the antennas is hopeless, we resorted to a sim¬ 
plified amplitude analysis that we have used before for 
fringe search surveys (Petrov at ah 2012a,b). Using the 
measured Tgyg (a modelled Tgyg in the case of WARK30M ) 
and the a priori gain curves, we computed a priori SEFD 
for each observation and calibrated fringe amplitudes av¬ 
eraged over time and frequency. In a case of CEDUNA only 
the fringe amplitude from the first IF (6.592-6.608 GHz) 
was used. We identified 48 sources with known images 
from the VLBI Image Database^ maintained by the As- 
trogeo Center. Using these images, we computed the 
expected fringe amplitude at given baseline projections. 
We assumed that the image is the same at 8.4 and 
6.7 GHz. We formed the ratios of the measured fringe 
amplitude to their expected values and solved for loga¬ 
rithms of antenna-based gain factors. Estimates of the 
SEED at 45° elevation angle are given in Table 2. The 
anomalously high SEED at HOBART is due to a failure of 
the cryogenic equipment that resulted in a system tem¬ 
perature rise to ^200 K during the experiment. 

TABLE 2 

Average fitted SEFD at 45° elevation angle. 

GEDUNA 890 Jy 

HOBART26 1890 Jy 
wark30m 650 Jy 

The error of SEFD estimates is 15%. 

Although the flux density of an individual source may 
change due to variability, the average ratio of the instan¬ 
taneous flux density at the experiment epoch to its value 
on the image epoch for the ensemble of 48 sources should 
be stable. Using a relatively large ensemble for the evalu¬ 
ation of gain factors allows us to estimate their statistical 
uncertainty: 15%. 


TABLE 3 

The median correlated flux density at 6.7 GHz 

DETERMINED FROM ANALYSIS OF THIS EXPERIMENT. 


lAU Name 

IVS Name 

(3) 

(4) 

(5) 

(6) 

J0004-4736 

0002-478 

c 

6 

0.38 

0.03 

J0012-3954 

0010-401 


3 

0.36 

0.03 

J0019+2021 

0017+200 

c 

3 

1.03 

0.04 

J0038-2459 

0035-252 


3 

0.71 

0.03 

J0049-5738 

0047-579 


3 

0.68 

0.03 

J0050-0929 

0048-097 

c 

3 

0.36 

0.03 

J0058-0539 

0055-059 

c 

3 

0.35 

0.02 

J0058-5659 

0056-572 


3 

0.44 

0.03 


(3): Flag “C” denotes the source was used as an amplitude calibra¬ 
tor; (4): The number of observations; (5): Correlated flux density 
in Jy; (6): uncertainty of correlated flux density in Jy. 

Table 3 is published in its entirety in the electronic edition of the 
journal. A portion is shown here for guidance regarding its form 
and content. 

® Available at http://astrogeo.org/vlbi_images 



Fig. 5.— Field around 1031-837 at 4.6 iim. from WISE. The 2- 
(7 error ellipse around 2FGL J1032.9-8401 position denoted as a 
triangle is shown with red colour. The blue dot shows 1031-837 
position. 

Using empirical gain factors, we evaluated the aver¬ 
age flux density of each of the 238 sources observed in 
this experiment. The average flux density of 1031-837 
was found to be 46 ± 7 mJy. Since it was detected only 
in four scans on the HOBART26/WARK30M baseline, no 
information about its structure can be derived. Never¬ 
theless, using information that we derived from analysis 
of this experiment, namely, its correlated flux density at 
projected baseline 50 MX and its position, 1.96 Fermi 
position uncertainty, we were able to compute the like¬ 
lihood ratio of its association with the y-ray object as 
defined in Schinzel et al. (2015): 22.8. That means that 
the probability of the y-ray source and the radio source 
being the same object is by a factor of 22.8 greater than 
the probability of two unrelated objects being inciden¬ 
tally projected close to each other. Thus, we conclude 
that we have established an association between the ra¬ 
dio and y-ray sources. 

There is a source WISE JI030I5.28-840308.6 from 
the Wide-held Infrared Survey Explorer (WISE) cat¬ 
alogue (ALLWISE, November 13, 2013 (Wright et al. 
2010; Mainzer et al. 2011)) within 90 mas of I03I-837 
position (see Eigure 5). Cross-matching the WISE and 
the Radio Eundamental Catalogue hnds 5984 matches. 
Eitting parameters of the Rayleigh distribution to dis¬ 
tances between source positions of matches we found that 
the estimate of the Icr uncertainty along each coordinate 
axis is 95 mas. We interpret this estimate as a measure of 
position accuracy of the ALLWISE catalogue for point¬ 
like sources. Therefore, we conclude that I03I-837 and 
JI030I5.28-840308.6 are the same object. It was detected 
in all 4 bands of WISE instrument with magnitudes of 
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14.10, 13.11, 10.27, and 8.16 at 3.4, 4.6, 12 and 22 /im, 
respectively. Figure 5 demonstrates that the number of 
IR sources within the Fermi position error ellipse is so 
high, that their association with the 7 -ray object is im¬ 
possible without use of additional information. 

Analyzing the flux densities of other sources, we found 
that the typical detection limit for 90 s integration on the 
HOBART26/WARK30M baseline ranged from 40 to 70 mJy 
with a median value of 50 mJy. Fixing the cryogenic 
problem experienced at HOBART26 promises to reduce 
the detection limit to a 20 mJy level. 

4. SUMMARY AND FUTURE WORK 

A new radio astronomy antenna WARK30M has ob¬ 
served successfully its first scientific VLBI experiment. 
Its position with uncertainty over the vertical coordinate 
5 cm and horizontal coordinates 1 cm has been deter¬ 
mined. Error analysis shows that the systematic errors 
are a factor of 2 smaller than these uncertainties. The 
averaged SEFD at 45° elevation angle was 650 Jy. 

The first scientific experiment allowed us to estab¬ 
lish an association of Fermi detected 7 -ray source 2 EGL 
J1032.9-8401 that was previously considered unassoci¬ 
ated. Position of its radio counterpart was determined 
with the uncertainties of 4.6 mas over both coordinate 
components and its correlated flux density was found to 
be 54 ± 5 mJy. We conclude that WARK30M is ready for 
scientific programs. 



Fig. 6. — The uv coverage of simulated 12 hour loug observatious 
of 1934-638 at a 7-statiou uetwork wark30m, ceduna, parkes, 
HOBART26, ATCA, MOPRA, aud HARTRAO. Poiuts at baseliues with 
WARk30m are showu with red. 

At the moment, there are several VLBI arrays that 
are actively observing: the Very Long Baseline Ar¬ 
ray (VLBA), the European VLBI Network (EVN), the 
VLBI Exploration for Radio Astronomy (VERA), the 
Korean VLBI Network (KVN), the Chinese VLBI Net¬ 
work (CVN), the Russian Quasar network, and the Long 
Baseline Array. All these arrays, except the last, are 
located in the northern hemisphere and cannot observe 


sources with declinations less than —40°. WARK30M will 
be a valuable addition to the network of sensitive south¬ 
ern telescopes (SEED less than 1000 Jy): 26m Hartrao, 
30m Ceduna, 64m Parkes, 22 m Mopra, 6 x 22 m ATCA, 
and 26 m Hobart. Eigure 6 shows the i^i^coverage when 
this seven-station array observes the 1934-638. Tracks 
of baselines with stations to WARK30M are highlighted 
in with red. The main focus of VLBI programs with 
WARK30M is anticipated to be participating in surveys 
targeting sources in the declination zone [—90°,—40°] 
and exploring peculiar objects, mainly in this declina¬ 
tion zone inaccessible to the arrays in the northern hemi¬ 
sphere. These programs will be complement to programs 
that are conducted in the northern hemisphere and will 
convert them into all-sky surveys. We anticipate that 
the most valuable contribution of the southern VLBI ar¬ 
ray with WARK30M participation will be given to those 
programs that benefit from completeness and all-sky cov¬ 
erage. These programs include the absolute astrometry 
survey that is used for space navigation, for a connec¬ 
tion of VLBI and Gaia coordinate systems, as a list of 
phase calibrators for phase referencing programs, and for 
other important applications, such as the VLBI survey 
of a flux-limited sample of compact radio sources; for a 
search for associations of 7 -ray sources. We have plans 
to equip WARK30M with a cooled C-band receiver, an X- 
band receiver and with the a 2 Cbps recording system in 
2015. 

Among the 244 observed sources, 175 have been de¬ 
tected in three or more observations. We determined 
their correlated flux densities with errors 15% (See Ta¬ 
ble 3). This shows the potential of a southern hemisphere 
VLBI array to contribute to programs of flux density 
monitoring of active galaxy nuclei associated with 7 -ray 
sources. The TANAMI program of monitoring 84 ra¬ 
dio sources associated with 7 -ray objects in the south¬ 
ern hemisphere with the Australian Long Baseline Ar¬ 
ray (LBA) (Ojha et al. 2010) whose goal is to provide 
a time series of images. However, TANAMI observes 
25-30 sources per 24 hour experiment at a 6-10 station 
network, i.e. requires more than one order of magnitude 
more resources. A program of coarse flux density mon¬ 
itoring of a much wider list will be complimentary to 
TANAMI. By March 2015, there are 270 known 7 -ray 
sources with declinations < —30° associated with the ra¬ 
dio sources that exhibit emission from milliarsec regions 
detected with VLBI. All these sources can be observed in 
a 27-30 hour observing session. A program of monthly 
monitoring of these sources has the potential to check 
their variability and derive light curves. 

Both VLBI and single-dish spectral observations in C- 
band can potentially include spectral line observations 
of OH maser lines at 6.03 and 6.035 CHz and methanol 
masers at 6.7 CHz, which are important for the iden¬ 
tification and study of star-forming regions. Eor exam¬ 
ple, 6.7 CHz methanol masers are exclusively associated 
with the early stages of high-mass star formation and 
VLBI observations can be used to obtain accurate dis¬ 
tance estimates through trigonometric parallax. Single 
dish observations with the 30m radio telescope in New 
Zealand can be used to monitoring of the variability of 
methanol maser sources found by the Parkes multi-beam 
survey (Woodburn et al. 2015). 
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